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^ : Abstract 

lO ! The exact energy spectrum is developed for a two temperature kinetic Ising spin chain, and 

its dual reaction diffusion system with spatially alternating pair annihilation and creation rates. 
' (— I '. 

Q ■ Symmetries of the system pseudo-Hamiltonian that enable calculation of the spectrum are also 

^ . 

^ , used to derive explicit state vectors for small system sizes, and to make observations regarding 



state vectors in the general case. Physical consequences of the surprisingly simple form for the 
■ eigenvalues are also discussed. 
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I. INTRODUCTION 



Over the last three decades, an increasing number of condensed matter theorists are 
devoting their efforts to understanding complex collective behavior of far-from-equilibrium 
systems using methods that range from easily accessible computer simulations to sophis- 
ticated theoretical studies. Although great progress has been made, a comprehensive 
theoretical framework is still lacking. In this context, low- dimensional systems are of 
particular interest since their simplicity permits analytical and numerical solutionsi'^i'^ 
and these non-trivial solutions shed light on related, more complicated higher-dimensional 
models. 

This paper presents the exact energy spectrum of two closely related one- dimensional 
non-equilibrium models: a kinetic Ising chain (KISC) with cells coupled alternately to one 
of two temperature baths, with generalized Glauber dynamics^, and its dual counterpart, a 
reaction diffusion system (RDS) with spatially alternating pair annihilation and creation 
rates. Interest in these models is also motivated by their experimental applications. 
Multi-temperature spin systems are fairly common: nuclear magnetic resonance in an 
external magnetic field is an example; a lattice of nuclei in a solid prepared at a finite spin 
temperature^ is another. On the other hand, the RDS model with spatially alternating 
annihilation and creation rates is known to describe the dynamics of photo-excited solitons 
in polymers^. Mobilia et al. proposed an experimental realization of the RDS model with 
alternating rates in MX chain compounds using a laser with spatially modulated power 
output^. 

The two-temperature kinetic Ising model (KISC), was first introduced by Racz and Zia^ 
who calculated exactly the two-point correlation functions for the steady state. Using a 
perturbation expansion of the master equation, Schmuser and Schmittmann^ calculated the 
first two corrections to the equilibrium Boltzmann distribution. Mobilia et al? found an 
analytical solution for the full dynamics (magnetization, particle density, and all correlation 
functions) of this non-equilibrium spin chain and its related reaction-diffusion model using a 
generating function approach. Outstanding challenges include knowledge of the exact energy 
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spectrum of these models, and a compact expression for their steady states. 

Our study brings us one step closer to achieving this goal. Using the standard mapping^ of 
reaction diffusion models onto integrable quantum chains, the RDS model can be expressed 
in a "free fermion" form, by defining a quadratic non-Hermitian "stochastic Hamiltonian"^. 
This operator can be diagonalized as long as certain constraints are obeyed^. In this paper, 
we derive the exact energy spectrum of this pseudo-Hamiltonian. We also utilize symmetry 
considerations to extract the energy eigenvalues and associated eigenstates for some small 
system sizes with the goal of finding a general pattern for the steady states of these models. 

Our paper is organized as follows: In Section 1, we give an overview of the models. Next 
(Section 2), we describe the symmetries exhibited by the pseudo-Hamiltonian operator, and 
their role in the diagonalization process. We present exact solutions (eigenvalues and eigen- 
vectors) for some small system sizes in Section 3. Following some standard technical steps 
(Jordan Wigner transformation, discrete Fourier transform and a generalized Bogoliubov 
transformation) we derive closed-form expressions for the eigenvalues, and a methodology 
for extracting the eigenvectors. Section 5 presents a summary of our results and some pos- 
sible generalizations for these models. 



II. OVERVIEW OF MODEL 



Two equivalent one-dimensional models motivate the work herein: the kinetic Ising spin 
chain (KISC) and its associated reaction-diffusion model (RDS.) The KISC model parallels 
the one-dimensional Ising model. We postulate a lattice of side-by-side cells, numbered 
n = 1,2, . . . , N, arranged in a ring such that cell n = is considered adjacent to cell 
n = 1. N is restricted to even values. Each cell has a single degree of freedom with two 
possible values: —1, that can be thought a cell occupied by a particle with spin down; and 
+1, describing a cell occupied by a particle with spin up. Each cell interacts with its two 
nearest neighbors, as well as being in contact with a heat bath at one of two temperatures 
— Tg for even-numbered cells and Tq for odd numbered cells. If Tg ^ To the system cannot 
achieve equilibrium: each heat bath tries to drive the system towards a different equilibrium 
state. As a result, energy flows continuously between the even cell sublattice and the odd. 
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Configuration C (a list of the states of the N cells) changes into a different configuration C 
with generalized Glauber transition rates r [C ^ C] . Rate r is non-zero only if C and C 
differ only in the spin of a single particle. The rate at which site n has its spin flipped is 
given by: 

rn^^ - '^dn{dn-l + C?n+l) (1) 

where ks is Boltzmann's constant, the factor 7^ (0 < 7^ < 1) is related to the temperature 
of cell n by 

tanh(^^r) n even 

7n = < for (2) 

tanh(^^^r) n odd, 

and dn is the state (+1 or —1) of the n-th cell. This rate equation prescribes a spin flip rate 
for a cell of 1/2 if cells to the left and right have opposite spins, (1 — 7n)/2 if adjacent spins 
are the same and the same as that of cell n, and (1 + 0n)/2 if adjacent spins are the same 
and opposite that of cell n. The time scale is arbitrary. 

The KISC model is mapped onto an equivalent reaction-diffusion model with spatially 
alternating pair creation and annihilation rates in the following way. A dual lattice of N 
sites is established, in which a site in the dual lattice is associted with the boundary between 
two sites in the KISC lattice. A pair of adjacent KISC spins with opposite signs is identified 
with a particle in the dual lattice; adjacent spins with the same sign is identified with the 
absence of a particle (a hole.) A spin flip in the KISC model translates into either diffusion 
of particles on the dual lattice with equal left-right rates, or pair creation or annihilation 
with different rates. Transition rates between configurations in the KISC system become 
diffusion, pair creation and annihilation rates in the RDS system, as shown in Table 1. 

Time evolution of these systems is described by the master equation, expressing conserva- 
tion of probability assuming a continuous-time dynamics. The probability P{C, t) of finding 
the system in configuration C at time t increases due to transfer of probability into C from 
other configurations, and decreases as C passes probability into others, in such a way that 
^^P(C, t) = 1 for all t. The evolution of probability P{C,t) is described by transition 
rates r [C — C] , the probability per unit time that configuration C changes into a different 
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configuration C . The master equation is: 

=Y,{r [C - C] P{C', t)-r[C^ C] P{C, t)} (3) 

C'^C 

in which the first term on the right represents the gain in probabihty of configuration C 
due to transitions from other configurations, and the second represents losses due to C 
transforming into other configurations. 

We utihze Dirac notation to represent each configuration as \C). From this we build 
a vector representation of a probabilistic superposition of all possible configurations of a 
system: 

\P{t)) = J2PiC,t)\C). (4) 

c 

The master equation can now be re-expressed in terms of this vector as: 

||P(t)) = -H\P{t)) (5) 

where the pseudo-Hamiltonian if is a 2^ x 2^ matrix, with matrix elements 

{C'\H\C) = -r{C ^C'),C' y^C (6) 
{C\H\C) = Y.r{C-^C'). (7) 

A formal solution to Eq. Ocan be written as \P{t)) = e~^'^\P{0)). Our goal is to investigate 
the eigenvalues of operator H in order to explain the system's time dependence. 

From this point, we shall focus on the RDS model. We follow the precedent of repre- 
senting the "particles" and "holes" in the dual lattice by a spin one-half model: a particle 
is represented by spin up (and thus a hole becomes a spin down (|0).) From the above 
formalism comes the definition of the probability-conserving operator H that controls the 
system's time dependence: 

- 2^ = 5^ + + (1 + 7e)^/^;+l + (1 - 7e)^7a7+l - 7e(^7a+ + aT^ia+^J 
jeven 

- (1 - 7e)] + + + (1 + 7o)fT+(T++i + (1 - 7o)^7^7+l 

jodd 

- lo{(r-al + - (1 - 7o)]. (8) 
The operators and cx^ are the Pauli spin raising and lowering operators on the n-th cell: 

a+|0)„ = cr+|l)„ = 0, a^|0)n = 0, and a^|l)n = |0)n. 



It has been shown that the eigenvalues and eigenvectors of H can be found if the "free 
fermion constraint" is obeyed^. For the RDS model, this means that the sum of local 
diffusion rates is equal to the sum of local pair creation and annihilation rates (in our case 
I + ^ = '^^2°'" ~^ ^^2°'" •) This constraint assures the biliniarity of the H operator, and, 
consequently, an exact solution for the problem. 

III. SYMMETRIES OF THE H OPERATOR 

Symmetries exhibited by the H operator, Eq. [8|, affect the form of its eigenvalues and 
eigenvectors, and in some cases aid in the process of determining them. In this section our 
goal is to exhibit these symmetries towards the goal of direct calculation of eigenvalues and 
eigenvectors of H for small values of A^. 

Because the fundamental process described by the rate equation, Eq. [H corresponds to 
the simultaneous flipping of two spins in the RDS model, the H operator does not change 
the 'spin-parity' of a state, i.e. states with an even number of up spins are transformed by H 
into states of only even number of up spins, and likewise for odd numbers of up spins. This 
symmetry immediately separates the configuration space into two sub-spaces of the same 
dimensionality (2^/2) which do not interact. Thus the 2^ x 2^ if matrix is reduced to two 
equal-sized diagonal blocks by proper ordering of the configuration basis states. 

The H matrix is also invariant to a translation of the ring of cells by an even number of cells 
to the right or left. Thus H commutes with the operator that invokes this translation, and 
simultaneous eigenstates of the two operators can be found. Such eigenstates are conveniently 
written as sums of the form 

N/2-1 
n=0 

{q = 0, . . . , N/2 — 1), where \Cq) is a spin configuration (e.g. |011000) representing six 
spins, the second and third ones up, for the 1x6 case) and |C„) is the same configuration 
pushed 2n cells to the right using periodic boundary conditions (thus if |Co) = |011000), 
then IC2) = 1 100001).) The translation symmetry of the pseudo-Hamiltonian implies that H 
does not mix states of the form \Cq) which have different q values. Thus within each of the 
two major sub-blocks of the H matrix (one of even spin parity, one of odd spin parity) there 
are N/2 smaller sub-blocks, each with a different value of q. This symmetry has enabled 
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relatively straightforward extraction of eigenvalues and eigenstates for even N values up to 
N — S (with a 256-dimensional configuration space.) 

Two additional symmetries are apparent from the form of H given above, that provide 
further information regarding the form of the eigenvalues and eigenvectors. 

The simultaneous translation of the spin chain by a single site (cell n becomes cell n + 1) 
along with the interchage of values ■je ^ lo leaves H invariant. If we use X to represent this 
transformation, it follows that X"^ = I. If Itjj) is an eigenstate of H with eigenvalue E, then 

This leads to several possibilities: 

• is invariant under the interchange je ^ lo and either X\%1)) is a constant multiple 
of \'4>), or X\'^) produces another eigenstate of H distinct from j'^) but with the same 
eigenvalue. 

• is not invariant under the interchange 7e <-> 7o, but instead transforms to another 
distinct eigenvalue of if, and X\%1)) becomes a corresponding eigenstate. 

Explicit diagonalization for small N , shown in the next section, suggests that the eigenvalues 
are invariant under 7e <-> 7o. The general solution for the eigenvalues presented in a later 
section show this to be true for all even values of N. 

Another symmetry operation leaving H invariant consists in changing the sign of both 7 
constants and simultaneously flipping every spin. If we use Xa to represent this transforma- 
tion, algebra similar to that of the prior paragraph leads to the following: 

• is invariant under the change in sign of both 7's and either Xjs\%l)) is a constant 
multiple of 1-0), or Xjs\%l)) produces another eigenstate of H distinct from l'^) but with 
the same eigenvalue. 

• is not invariant under the change in sign of both 7's, but instead transforms to 
another distinct eigenvalue of and Xji\%l)) becomes a corresponding eigenstate. 

Explicit diagonahzation for small N (following section) suggests also that the eigenvalues 
are invariant under the change in sign of both 7's. This, too, is later proven true for all N . 
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IV. EXPLICIT DIAGONALIZATION OF H 



The number of configurations for a given value of N ^ thus the dimensionahty of the square 
matrix i?, is 2-^. This prohibits exphcit diagonahzation of H for all but the smallest values 

of N . The symmetries of H enable straightforward use of computer algebra software to 
find eigenvectors and eigenvalues for several cases. In this section some explicit results for 
A?" = 2, 4, 6, 8, so calculated, are exhibited and discussed. 

A. N = 2 case 

The simplest case, N — 2, while simple is not quite trivial, and its solutions are useful 
for establishing patterns for eigenvalues and eigenvectors of H. For this case, there are four 
basic configurations (|00), |11), |01), |10)) and thus the H matrix is 4 x 4. The operation 
of shifting each configuration two cells to the right is equivalent to the identity operation, 
and this plays no useful role for this case. Spin parity does play a role, indicating that H 
simplifies into two 2x2 diagonal blocks if the basis configurations are ordered as above (even 
spin-parity states |00) and |11) form a basis for the even subspace, odd spin-parity states 
1 01) and 1 10) form a basis for the odd subspace.) The H matrix for this basis is: 



/ 1 I 7e+To 



_l + Ie±To 



\ 



H = 



-1 - 



7e+7o 

2 



1 



7e+7o 

2 



















1 



1 




-11/ 
0, the ground state with eigenvector 



(2-7e-7o){00} + (2 + 7, + 7„){ll}; 



and eigenvalue 2 with eigenvector 



{00} -{11}. 



The lower block also has two: eigenvalue with eigenvector 



{01} + {10}; 



8 



and eigenvalue 2 with normalized eigenvector 

{01} -{10}. 

The construction of H guarantees that its eigenvalues are non-negative, and that its lowest 
eigenvalue is zero. These results show a twofold degeneracy of the zero eigenvalue, with one 
eigenstate in the even subspace, and the other in the odd subspace. This happens as well in 
the following special cases presented below, and is a general property for all N . 

B. iV = 4 case 

This case requires the diagonalization of a 16 x 16 matrix, so utilization of symmetries of 
H is useful, if not essential. Both the spin parity symmetry and the translation symmetry 
play a role in this case. H breaks into two 8x8 diagonal blocks immediately - one of even 
spin parity states, the other of odd. Within each block, there are two smaller blocks when 
the basis states are transformed as in Eq. [9l We characterize basis states in the notation of 
Eq. [H] by their g- values. Of the eight even spin-parity states, six are g = states, 

^ |0000) \ 
|1111) 
|0101) 
|1010) 

|1100) + looii) 

\ |1001) + lOllO) y 
and the remaining two are g = 1 states, 

1001) - |0110) \ 
1100) - |0011) j ' 

Each of these separately are bases for a diagonal block of if, one 6x6, another 2x2. 

The 8x8 even-parity subspace block for H, using these eight basis states, g = states 
first, is 
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V 



2 + 7e + 7o 00 

2-7e-7<, 

2 

2 

-l-7e -l+7e -1 -1 

-I --to -l+7o -1 -1 







-l+7e 



-l+7o 







2 2 

-1^ ^ 



2 








1 










2 
2 



The upper 6x6 block (g = 0) of this matrix has eigenvalues 0,2,2,4,2 — y^27^, and 
2 + -\/27e7o- Corresponding eigenvectors are displayed in Fig. 2. 

Of the remaining eight odd spin-parity basis states, four can be cast as g = states, with 
basis vectors 

( |0010) + |1000)\ 
|1101) + loiii) 

lOOOl) + lOlOO) 
\ |1011) + |1110) ) 
and four as g = 1 states, with basis vectors 



/ lOOlO) - 


|1000) \ 


|1101) - 


lOlll) 


|0001) - 


|0100) 


Uoii)- 


11110) J 



Using this basis, with the g = states in the first four positions, this 8x8 diagonal block 
of the Hamiltonian from the odd spin parity basis breaks into a 4 x 4 block and two 2x2 
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blocks: 



2 + 





-1 


1 1 7e+7o 

2 

















O 7e+7o 

2 


1 7e +7o 

2 


-1 














-1 


-1 


2 + 

















1 7e+7o 

2 


-1 





O 7e+7o 

2 


























2 + ^^+^° 


-7e+7o 

2 




















7e-7o 

2 


O 7e+7o 

^ 2 


























2 + 


7e-7o 

2 




















-7e+7o 

2 


O 7.! +7 
^ 2 



The 4x4 block has eigenvalues 0, 2, 2, and 4; each of the succeeding 2x2 blocks has 
eigenvalues 2 ± ^^elo- Corresponding eigenvectors are displayed in Fig. 2. 

Note that as was the case for N = 2, there is a twofold degenerate ground state (eigenvalue 
0) with one eigenstate in the even subspace and the other in the odd subspace. 

C. A'^ = 6 case 

The configuration space for this case is of dimension 2^ = 64. As suggested by the results 
of the = 4 case, the form of the eigenvalues is quite complicated, and will not be displayed 
herein. The following discussion outlines the steps in solving for eigenvalues and eigenvectors, 
and exhibits the eigenvalues. 

H breaks into two 32 x 32 diagonal blocks - one of even spin-parity states, the other of odd. 
Within each block, there are now three smaller blocks when the basis states are replaced by 
linear combinations as in Eq. [9l Eigenstates in the respective sub-block bases are of the forms 
\ijklmn) + \klmnij) + \mnijkl) (sub-block 1, q = 0), \ijklmn) + Iklmnij) + e^~\mniikl) 
(sub-block 2, q = 1), and \ijklmn) + \klmnij) + \mnij kl) (sub-block 3, g = 2.) Note 
that there are four states, namely |000000), |010101), |101010), and |111111) that transform 
into themselves under forward translation by two cells, and thus "interact" (through H) only 
with other states in the q = sub-block. The first and fourth of these states are of even spin 
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parity and are part of a 12 x 12 block of even parity states with g = 0: 



000000) 

mill) 



110000) 


+ 


001100) 


+ 


000011) 


011000) 


+ 


100001) 


+ 


011000) 


001111) 


+ 


110011) 


+ 


111100) 


100111) 


+ 


111001) 


+ 


011110) 


001001) 


+ 


010010) 


+ 


100100) 


010001) 


+ 


010100) 


+ 


000101) 


101000) 


+ 


001010) 


+ 


100010) 


101110) 


+ 


101011) 


+ 


111010) 


110101) 


+ 


011101) 


+ 


010111) 


110110) 


+ 


101101) 


+ 


011011) 



Commercially available algebraic software was used to diagonalize this and the other 
blocks for this case (as for the previous cases.) The twelve eigenvalues for this block are 

E = (the even spin parity ground state), E = 2 with degeneracy 3, E = 4 with degeneracy 
3, E = 6, E = 2± V^^o, and £; = 4 ± ^StTt^ 

The remaining two states that transfer into themselves under translation by two sites — 
(|010101), |101010)) — are part of the q — sub-block for odd spin-parity. The twelve states 
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loioioi) ' 

llOlOlO) 
loioooo) + lOOOlOO) + lOOOOOl) 
liooooo) + looiooo) + lOOOOlO) 
liooioi) + lOllOOl) + loioiio) 

101001) + lOlllOlO) + IIOOIIO) 
IllOOOl) + lOlllOO) + lOOOlll) 

IllOlOO) + looiioi) + loiooii) 
IllOOlO) + llOllOO) + looioii) 
111000) + looiiio) + llOOOll), 

IlllllO) + llOllll) + IlllOll) 

IllllOl) + lOlllll) + IllOlll) y 

and the eigenvalues for this block are E — (the odd spin parity ground state) , E — 2 with 
degeneracy E — A with degeneracy 3, E — 6, E — 2± -sJ^^TTo-, ^'^^ -B = 4 ± ^/tT)^- 

There are two remaining even spin-parity sub-blocks, both 10 x 10 and both have a basis 
set of the form 



/ 


110000) + n 


001100) +n2 


000011) ^ 




011000) + n 


100001) 


011000) 




001111) + n 


110011) 


111100) 




100111) +n 


111001) +71^ 


011110) 




001001) +n 


010010) + 


100100) 




010001) +n 


010010) 


100100) 




101000) +n 


001010) 


100010) 




101110) +n 


101011) 


111010) 




110101) + n 


011101) 


010111) 


\ 


110110) +n 


101101) +71^ 


011011) y 



This is the even spin-parity sub-block for g = 1 when n — e*^'^/^, and even spin-parity sub- 
block for q = 2 when n = e*^'^/'^. For either of these blocks, the family of eigenvalues is the 
same: E = 2, E = 4:, E = 2± -\/37e7o/4 each with degeneracy 2, and E = A± -\/37e7o/4 
each with degeneracy 2. 



are 



v 
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The final two blocks are odd parity blocks, each 10 x 10, with basis sets of the form: 



/ 


n 1 r\r\r\r\\ i ^ 
UlUUUU/ + n 


UUUIUU/ + n 


UUUUUl; 




1 r\r\r\r\r\\ i 

100000) + n 


001000) + n 


000010) 




100101) + n 


011001) + n 


A1 A1 1 a\ 

010110) 




101001; + n 


011010) + n 


100110) 




110001) + n 


011100) + n 


A A A 1 1 1 \ 

000111) 




110100) + n 


001101) + n 


A i A A 1 1 \ 

010011) 




110010) +n 


101100) 


001011) 




111000) +n 


001110) 


100011) 




111110) +n 


101111) 


111011) 


V 


111101) +n 


011111) +71^ 


110111) J 



This is the odd spin parity sub-block for g = 1 when n — e*^^/^, and odd spin parity sub- 
block for = 2 when n — e'^'^/^. For either of these blocks, the family of eigenvalues is 
the same: E^2, E^A, E^2± E ^ A± ^/77)^/2 , E ^2± 3^/77)^/2 and 

E — Ark 3i/7e7o/2 each with degeneracy 2. 

In a similar, but much more tedious procedure, the eigenvalues and eigenvectors of the 
N — 8, 256 X 256 H matrix can be found explicitly The results for that case and the 
smaller- A?" cases support the following claims, to be proven for general N in the following 
section: 

• there are exactly two eigenvalue states; one in the even and one in the odd spin-parity 
subspace; 

• the maximum eigenvalue in both the even and odd spin-parity blocks has value N; 

• each eigenvalue state is in the q — subspace, with eigenvectors containing symmetric 
sums of shifted states; 

• eigenvalues depend only upon the single parameter ^J^y7To■ This implies that the 
spectrum can be fully deduced through knowledge of the spectrum for the single- 
temperature case, Tg = To- 
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V. GENERAL CASE 



Here we approach the problem using standard methodology^. Starting with the full, 
two-temperature operator in terms of spin raising and lowering operators, Eq. [8], we can 
rewrite in terms of fermionic operators {cj, cj) that satisfy anti-commutation relations 



{cj,cl} = Cjcl + clcj = 6. 



{cj,Ci} = {c],cj} = 0. 
This transformation, due to Jordan and Wigner^^ is as follows 



0" = Cje~*''^"<-' 



Straightforward application of this transformation produces the 'fermionized' pseudo- 
Hamiltonian 



j even 

- (1 + 7e)] -\Y1 + SVlCj + (1 + 7o)c]c]+i - (1 - -fo)CjCj+i 

j odd 

+ ^o{c]c, + cj+ic,+i) - (1 + 7o)]. (10) 



For the even spin-parity subspace the requirement that cr^^i = cr^ and cr^^^ = cr]" 
implies cn+i = — ci and cjv+;^ = —c\. For the odd spin-parity case, cn+i = ci and cjy+i = c\. 
We take advantage of the translation symmetry mentioned in Section 3 by defining two 
kinds of fermions in momentum space, one created from even-numbered cells, the other from 
odd-numbered cells. Define the following momentum-space operators: 



m even Q&Q 
j odd g&Q 

We chose q values to belong to the set Q = {±^, . . . ± -^^^^r^} for states with 

even spin-parity, and Q = {0, ±7^, i^f , . . . ± -^^^^p^, tt} for states with odd spin-parity, to 
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assure proper periodic boundary conditions for each case. These definitions of Q assume 
that N/A is integer- valued, but the uhimate results are valid as long as is even. 

Like the operators Cj and cj, the momentum space operators obey the canonical fermionic 
anticommutation relationships: 

{aq,aq>} = {4,4,} = {aq,al>} = ^q,q' 

= {6j,6;,} = {hq,h\,} = 5,,,, 

{a,,V} = K,6;,} = {S,&i,} = 

In terms of these operators, the pseudo-Hamiltonian is written as 

H = Y, \-coA{e''^a\hq + e-'h\aq) 



'2' 

+ alblq ^[(1 + le)e'' - (1 + 7o)] + a,6_, ^[(1 - 7e)e-^^ - (1 - 7o)] 
- I (4S + bib,) + 2 + 7e + 7o 



(11) 



The ultimate step in the derivation is to perform a Bogoliubov type similarity transform 
to new variables in which H takes diagonal form. We postulate a diagonal form for H that 
reads 

H = ^{uJqXqXq + ^'q^q^q + COUSt.), (12) 

while assuming that the operators Xq ^ind C,q obey Fermionic anticommutation relations. 
Since H is not Hermitian, Xq xl^ iqi^Ci- easily demonstrated that 

[Xq, H]^= UJqXq aud [^g, H]_= Uj'^^q. (13) 

Following Lieb, Schultz and Mattis^ we define the Xq and variables in terms of the a, 
and bq variables in the following manner: 

Xq = cial + C2b-q + c^ib^ + C4a_g (14) 

and 

^q = dial + d2b^q + d^bl + d^a^g, (15) 

where the q's and (ij's are constants. We can calculate the commutator with H in the form 
given by Eq. [TT] of each term in Xq as expressed in Eq. [HI using 

[at, H]^ = 7,„at + e'6_, + cos |e-^^/2 
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-lava- 



lav b\ 



Constants in these expressions are defined as: 

7e + 7o 



+ ebi - cos ?e-^«/2 ^_ 



e'a_o + cos ^e^"^'^ al 



t ^ /a/2 

ea^ — cos -e ^' a_q. 



la 



Idif 



,iq/2 



+ lav) sm - + iidif cos 



le-lo 

2 

-iq/2 I Q 
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e represents e as defined above but with q replaced by — g, and hkewise for e'. Eq. [13], with 
these results used to evaluate the left hand side, produces a large operator equation. The 
coefficients of each of the four operators on the right-hand-side must equal the corresponding 
coefficient on the left-hand-side. This requirement produces a linear set of equations for the 
constants q that can be expressed by the matrix equation UqC = MC where 
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Obviously, the eigenvalues of M are possible values of the energies {uq) of the x excitations, 
which is what we seek. Identical analysis of the di coefficients show that the eigenvalues of 
the same matrix are possible values of the energies {uj'^ of the ^ excitations. The eigenvector 
components are quite complicated, but the eigenvalues simplify easily. Two of the eigenvalue 
are negative, thus inappropriate candidates for the excitation energies (that we know to be 
positive.) The other two eigenvalues are 1 ± cos ^^Jlelo■ Identify a;^ = 1 + cos ^^Jl7io and 

= 1 — cos |^7e7o in order to deduce the exact spectrum for the model. 

There is a unique even spin-parity vacuum state |0)e defined by the relations: 

2)7r. 



X,|0)e = 0, e,|0)e = for gG{±^,±^, 



lOvr (A^ 



N 



■}• 
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Other even spin-parity states are formed by an even number of excitations of the Xq '^^ 
type: each excitation of the former type carries energy Uq, and each of the latter type carries 
ujg. Since these excitations are fermionic, there cannot be two Xq excitations with the same 
q, nor can there be two excitations with the same q. The highest energy state has N/2 
distinct Xq excitations and N/2 distinct C,q excitations, and carries total energy N. 

The odd spin-parity sector of the spectrum also has a unique vacuum: 

477 Svr (N-A)^ 
Xq\0)o = 0, ^g\0)o = for ; g e {0, ±— , ±— , . . . ± , 

Other odd-parity states are formed by an even number of excitations of the Xq or ^q type 
with q with values from the list appropriate to odd spin-parity. The highest energy odd 
spin-parity state has N/2 distinct Xq excitations and N/2 distinct C,q excitations, and carries 
total energy A^. 



VI. CONSEQUENCES AND CONCLUSIONS 



It is surprising that a relatively simple set of eigenvalues emerges from the great algebraic 
complexity of the solution for arbitrary of the two-temperature model considered herein. 
The eigenvectors exhibited in Fig. 2 for the = 4 case show that even in this relatively 
simple case, the state vectors are algebraically difficult. In particular, the E = eigenstate 
that corresponds to the steady-state solution, is not easily characterized. In general the 
eigenstates depend separately upon the values of 7e and 7o. While the methodology employed 
in the previous section can, in principle, allow the extraction of the eigenvectors of H, the 
algebraic complexity of the eigenvectors of the matrix M suggest that such a straightforward 
exposition of them is not likely to be illuminating. Nonetheless, our results produce numerical 
expressions for the eigenstates given specific values for the two temperatures, and thus to 
expressions for particle densities and correlation functions. General relationships for these 
have been previously exhibited by Mobilia et al^ 

The fact that the eigenvalues depend upon the single parameter ^J'^elo allows some simple 
deductions regarding special cases of the two-temperature model. If one of the temperature 
baths has infinite temperature (e.g. 7e = 0), the eigenvalues are those of the case of a single- 
temperature model with both baths at infinite temperature. These eigenvalue are identical 
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to those of the Glauber model^, but the eigenvectors are of greater complexity. If one of the 
temperature baths has T = (e.g. 7e = 1 ), the energy spectrum becomes the same as the for 
the one-temperature case related to the temperature of the other bath. In the RDS language, 
this case corresponds to a system with pair creation prohibited, and pair annihilation at a 
rate of 1 for the even sites. In general, for every case with distinct temperatures Tg and T^, 
there is a single temperature that will yield the same energy eigenvalues. The spectrum of 
energies for the single temperature case follows from the work of Grynberg et al}- (for the 
special case h = h'= i, e = (1 + 7)/2, e'= (1 — 7)/2 with 7 = tanh (2/A;bT). 

A few general observations about the spectrum in the thermodynamic limit, N ^ 00, 
are possible. Paralleling an observation of Grynberg, et. al.- for the single temperature 
case, as long as ^7e7o < 1 (at least one temperature bath is above absolute zero), there 
is a gap between the ground state and the next-highest energy level of 2(1 — ^/yTfo)- This 
assures that states other than the steady state decay exponentially in time. The spectrum 
of remaining states consist of bands of energy levels centering on states with E = 4,6,8, .. .. 
The widths of these bands grows with E, while the spacing between adjacent states remains 
constant at 2. As a result, regardless of how small the parameter y/jelo is, the bands will 
overlap for high energies. These observations are separately true for the even spin-parity 
and odd spin-parity segments of the energy spectrum. 

Because H does not cause transitions between states of different spin-parity, the time 
evolution of any initial state can be broken into two independent segments. Any initial 
configuration of spin states can be broken into a piece with even spin-parity with probability 
Pe and one of odd spin-parity with probability Po, with P^ + Po = 1. The even spin-parity 
segment decays towards the steady state |0)e maintaining constant probability Pe] likewise 
the odd segment decays towards |0)o maintaining constant probability Pq. The separate even 
and odd spin-parity energies control the rate of decay of the non-steady-state components 
for each segment. 

Symmetries discussed in Section 2 enable a few additional comments about the general 
form of the eigenvectors. Because we know the eigenvalues to be invariant under the in- 
terchange 7e ^ 7o it follows that displacement of the ring by a single site (effectively a 
permutation of basis states) along with an interchange of values of 7e and 70 should trans- 
form any eigenvector into a constant multiple of itself, or into a different eigenvector with the 
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same eigenvalue. Since the ground state {E = 0) and the maximum energy state {E = N) 
for the even spin-parity sector and for the odd spin-parity sector are nondegenerate, each 
should be invariant within a constant under this transformation. 

The eigenvalues are also invariant under the simultaneous change of sign of 7e and 7o. This 
implies that a simultaneous flip of all spins (another basis state permutation) accompanied 
by a sign change in both 7s should transform an eigenvector into a constant multiple of itself 
or into another eigenvector with the same E. The non-degenerate eigenstates in each spin- 
parity sector should transform into constant multiples of themselves under this symmetry 
transformation. 

Mobilia et air have examined the behavior of this model in the case where the 7s have 
opposite signs. Although the concept of negative temperatures does not make physical sense 
for the KISC model, in the context of the RDS system it corresponds to a grid where the 
pair creation rate exceeds the pair annihilation rate on one sub-lattice, and the opposite 
is true for the other sub-lattice. The eigenvalues for this case will have positive real parts 
{Re{E) = 0, 2, 4, ... , N) corresponding to exponential damping in time, and imaginary parts 
proportional to A/|7e7o| producing oscillatory behavior. This result is consistent with the 
predictions presented in Mobilia et al.^, for example that under these conditions the density of 
particles approaches its equilibrium value via a term proportional to exp(— 2t) sm{2^/[y^rf^t+ 
6). Our spectrum results show such a behavior for the highest frequency oscillation associated 
with the most-slowly decaying component of the state function evolving from a general initial 
condition. 

In ongoing work we seek a compact expression for the steady state of these models. Given 
the relationship between spin systems and reaction-diffusion systems, it will be interesting 
to investigate the effect of various initial conditions and open boundary conditions on the 
dynamics of the system. From an experimental point of view, open boundary conditions for 
RDS systems would be important in the study of chemical reactions that include creation 
and annihilation processes and dimer deposition. Although particle densities and correlation 
functions can be calculated fairly straightforward in the thermodynamic limit, the finite-size 
effects may also be worth investigating. 

We can also imagine other extensions of the models presented. For example, we are 
interested in considering an RDS model with non-uniform diffusion rates for the odd and 
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even sites, and different creation and annihilation rates. This can also shed some light on 
the general problem of dimerized spin chains^^. 
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FIG. 1: Correspondence between the two temperature kinetic Ising chain (KISC) and the equivalent 
reaction-diffusion model (RDS.) 
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Ordered basis states for even spin-parity: 
I0000>, lllll>, I0101>, I1010>, I1100>, I0011>, I1001>, I0110> 



even spin-parity eigenvalues, eigenvectors. N = 4 
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Ordered basis states for odd spin-parity: 
I0010>, I1000>, I1101>, I0111>, I0001>, I0100>, I1011>, Illi0> 



odd spin-parity eigenvalues, eigenvectors, N = 4 
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FIG. 2: Eigenstates (un-normalized) for the N = 4 case of the two-temperature kinetic Ising model. 
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